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Experimenters in Acoustics have discovered more than one set of phenomena 
apparently depending for their explanation upon the existence of regular currents of 
air, resulting from vibratory motion, of which theory has as yet rendered no account. 
This is not, perhaps, a matter for surprise, when we consider that such currents, 
involving as they do circulation of the fluid, could not arise in the absence of friction, 
however great the extent of vibration. And even when we are prepared to include in 
our investigations the influence of friction, by which the motion of fluid in the neigh- 
bourhood of solid bodies may be greatly modified, we have no chance of reaching an 
explanation, if, as is usual, we limit ourselves to the supposition of infinitely small 
motion and neglect the squares and higher powers of the mathematical symbols by 
which it is expressed. 

In the present paper three problems of this kind are considered, two of which are 
illustrative of phenomena observed by Faraday.* In these problems the fluid may 
be treated as incompressible. The more important of them relates to the currents 
generated over a vibrating plate, arranged as in Chladni's experiments. It was 
discovered by Savart that very fine powder does not collect itself at the nodal lines, 
as does sand in the production of Chladni's figures, but gathers itself into a cloud 
which, after hovering for a time, settles itself over the places of maximum vibration. 
This was traced by Faraday to the action of currents of air, rising from the plate at 

* " On a Peculiar Class of Acoustical Figures; and on certain Forms assumed by groups of particles 
upon Vibrating Elastic Surfaces,'' Phil. Trans., 1831, p. 299. 
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2 LORD RAYLEIGH ON THE CIRCULATION 

the places of maximum vibration, and falling back to it at the nodes. In a vacuum 
the phenomena observed by Savart do not take place, all kinds of powder collecting 
at the nodes. In the investigation of this, as of the other problems, the motion is 
supposed to take place in two dimensions. 

It is probable that the colour phenomena observed by Sedley Taylor* on liquid 
films under the action of sonorous vibrations are to be referred to the operation of the 
aerial vortices here investigated. In a memoir on the colours of the soap-bubble,t 
Brewster has described the peculiar arrangements of colour accompanied by whirling 
motions, caused by the impact of a gentle current of air. In Mr. Taylor's experiments 
the film probably divides itself into vibrating sections, associated with which will be 
aerial vortices reacting laterally upon the film. 

The third problem relates to the air currents observed by Dvorak in a Kundt's 
tube, to which is apparently due the formation of the dust figures. In this case we 
are obliged to take into account the compressibility of the fluid. 

[My best thanks are due to Mr. W. M. Hicks, who has been good enough to 
examine the mathematical work of the paper. The results are thus put forward with 
greater confidence than I could otherwise have felt.] 

§ 1. In the usual notation the equations of motion in two dimensions are 
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and since the fluid is incompressible, 
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* Proc, Roy. Soc, 1878. 

f Edinburgh Transactions, 1866-67. 
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and 

du dv . 

ay ax 
so that 

vV _Il vV= !^v 2 f + *lv 3 ± (4)> 

v dt v dx v dy 

For the first approximation we neglect the right-hand member of (4), as being of 
the second order in the velocities, and take simply 

^{^-li^= < 5 )- 

The solution of (5) may be written" 55 * 

^ == ^] i 4"^2 (6), 



where 

1 d 
v dt 



V 2 ^=0, v 2 --^W a =0 (7). 



We will now introduce the suppositions that the motion is periodic with respect to 
x, and also (to a first approximation) with respect to t. We thus assume that \jj x and 
i/f a are proportional to cos kx, and also to e in K The wave-length (X) along x is 27r/k, 
and the period r is 27r/n. The equations (7) now become 

'|-* 8 )*=°. (|.-* 8 -!>*=° o». 

by which i/^ and \// 2 are to be determined as functions of y. If we write 



fc — fc —j— \ J? 



we have as the most general solutions of (8) 

t/, 1= A<r^+Be + ^ (10), 

i/r 3 =Ce-^+De + ^ (1]). 

With respect to the value of &', we see from (9) that it is complex. If we write 

P=P 2 cos 2ot 3 -=P a sin 2a 5 

then 

k'=~P cos &+? P sin a. 

* Stokes " On Pendulums," Camb. Phil. Trans., vol. ix., 1850, 
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4 LORD RAYLEICH ON THE CIRCULATION 

In all the applications that we shall have occasion to make, an approximate value 
of k' is admissible. On account of the smallness of v, n/v is very large in comparison 
with P, that is to say, the thickness of the stratum through which the tangential 
motion can be propagated in time r is very small relatively to the wave-length X. 
We may therefore neglect & 4 in the equation 

and take simply 

Again 

(sin a— cos a) 3 =l— sin 2a=^— -, 

so that the difference between cos a and sin a may be neglected. We will therefore 
write 

k'=/3(l+i) . . . .• (12), 

where 

0=v£' (13). 

We must now distinguish the cases which we have to investigate. In the first we 
suppose that a wave motion is in progress in a vessel whose horizontal bottom 
occupies a fixed plane y=0. We may conceive the fluid to be water vibrating in 
stationary waves under the action of gravity, the question being to examine the 
influence of the bottom upon the motion. If there are no other solids in the 
neighbourhood of the bottom, we may put D = 0, y being measured upwards, and 
/3 being taken positive. 

The conditions to be satisfied at y=0 are that u and v should there vanish. Thus 

A+B+C=0, -&A+&B-&'C=0, 

so that 



\jj=C\ — cosh faj+T sinh ky-\-e ^ L 



and 

u=G{ — h sinh ky-{-k' cosh Jcy—k'e"^}. 

At a short distance from the bottom, u=k'C. If we denote by u Q the maximum 
value of u near the bottom, we have 

k'G=u e ini cosJcx, 

and then 

. . 7 f coshh/ , sinhfa/ t e~ hy ] , N 

\ff=u e m cos kx< — , ' + — b~^~*~~17~\ • ••••• (14), 
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h . 



u=u e int cos hc\ — , rsinh %+ cosh Jcy—e l ' y \ (15), 



k 



v=u e int sin Jbs-i — p cosh %+ s ^ n -h %+ v 6 ""*^ 



. . (16). 



These are the symbolical values. If we throw away the imaginary parts, we have 
as the solution in real quantities by (12), 



\jj=u cos hx 



cosh hy 



sinii hv c 1 

cos (nt— £ir)H ^ cos ^*+^— ^ cos (/&«— ^tt— 0y) J . (17), 



U = U Q COS JcX 



v=UnSm Jcx< 



'0 



h sinh hy 



7 cos (n£ — J w) + cosh hy cos n< — e"^ cos (nt — /&/) J- . . (18), 



a J' cos (»<— £ir)+ siah Ay cos n t+j~^ cos (w£— iir— #/) 



(19). 



This is the solution to a first approximation. At a very small distance from the 
bottom the terms in e~^ become insensible. 

Although the values of u and v in (18) and (19) are strictly periodic, it is proper 
to notice that the same property does not attach to the motions thereby defined of 
the particles of the fluid. In our notation u is not the velocity of any particular 
particle of the fluid, but of the particle, whichever it may be, that at the moment 
under consideration occupies the point x, y. If sc+£ y-\-r) be the actual position at 
time t of the particle whose mean position during several vibrations is x, y, then the 
real velocities of the particle at time t are not u y v, but 

du j. du dv j. dv 

1 dor dy 1 l chr dy' ' 



and thus the mean velocity parallel to x is not necessarily zero, but is equal to the 
mean value of 



in which again 



du £ du 
dx dy '* 



g= u dt, y)=\vdt. 



f (10 1 

From the general form of u, viz., cos kx F(^/, t) } it follows readily that — £dt=0. 

For the second term we must calculate from the actual values as given in (18), (19). 
Thus 
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u sin kx\ h cosh ky . , , , , . . , , . k e~^ J . , . _ , 

V= — ~^~~\ J^T sm (^ — i^)+ smh % sm ^+^72 sin {nt—]pr—Py) 



dn 7 

— =t& n cos feci 



W cosh % 



^ cos (nt -— ^77) + A; sinh ky cos nt-{-y/2.f3e Py cos (n£ + ^7r — /3z/) >• , 



of which, the two first terms may be neglected relatively to the third (containing the 

dt/ 

large factor /3). The product of rj and — will consist of two parts, the first indepen- 

ay 

dent of t, and the second harmonic functions of 2nt. It is with the first only that we 
are here concerned. The mean value of the velocity parallel to x is thus 



u { ? sin 2hx e Py 
4t% 



k cosh hy cos fiy-j- x /2.fi sinh hy sin {JPy—\Tr)—ke~ Py 



On account of the factor e fSy , this quantity is insensible except when hy is extremely 
small. We may therefore write it 

/ii 2 cnv} 2n"V P \ 

-* ~ loos fiy+Py (sin Py— cos/&/) — e~ Pv V ■ .... (20), 

V (equal to kjn) being the velocity of propagation of waves corresponding to k and n. 

The only approximation employed in the derivation of (15) and (16) is the neglect 
of the right hand member of (4), and the corresponding real values of u and v could if 
necessary be readily exhibited without the use of a merely approximate value of k'. 
To proceed further we must calculate the value of 

n dy^ v d^f 

v ax v ay v ' 

in (4), for which it will be sufficient to take the values given by the first approxima- 
tion. Thus 

and by (17) 

d-\fr 2 nu cos kx e~ Py . , 

from which we find as the value of (21), 



nkUf? sin 2kx e Py 



'0 



4:v*/3x/2 



—-— — j sinh hy sin fiy — ^/2 cosh ky cos fiy-\- ^/2 e~ Py 

+ terms in 2nt. 
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On account of the factor e~ Py this quantity is sensible only when y is very small. 
We may write it with sufficient approximation 



nk%i^ sin 27cx e ^ 



—Pysbifiy—GOBfiy+e-toX (22). 



The terms in 2nt, corresponding to motions of half the original period, are not 
required for our purpose, which is to investigate the non-periodic motion of the second 
order. The equation with which we have to proceed is found by equating (22) to v 4 ^- 
The solution will consist of two parts, one resulting from the direct integration of (22) 
and involving the factor e" Pv , the second a complementary function with arbitrary 
coefficients satisfying v 4 */>=0. In the calculation of the first part we may identify v 4 
with d^/dy 4 *, on account of the smallness of h relatively to /3. In this way our equation 
becomes 

d^yjr nkiir? sin 2kx e~^ y f _ . _ _ , . 1 , ■ 

d(M*~ w&~~ L ~& y sin & y ~ cos Pv+ e \ • • • • ( 23 )> 

of which the solution is 






fcos^+isin^y+^ysin^+^e-^l . . (24). 



The complementary function, being proportional to sin 2kx, may be written 

^^ {(A+ B 3/ )e-^+(^+BV)e^}. 

If the fluid be uninterrupted by a free surface, or otherwise, within distances for 
which hy is sensible, we must suppose (A / +B / 2/) = 0, so that by (13) the complementaiy 
function may be written 

The condition that v (equal to — dxjj/dx) must vanish when y=0, gives A= — yf. 
For the velocity parallel to x we have 



u=^ 



U(? sin 2kx 



V 



"e &{ — sin /&/— \ cos fiy+±fiy cos fly— % fly sin fty—ie^} 

+/3- 1 e-^{B-~2Jc(A+By)}~]. 



In order that ^ should vanish when ;*/=0, w T e must have 

B=2AA+fi8=ti8-^A = fi8, 

approximately. Thus 
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u= JL —^ [f f %{ — sin $y—\ cos $y+\$y cos $y—\$y sin fiy—ie~^} 

+ | 6 -%{l-2%}] . (25), 
and 

2k ^p^\e~^{l cos #/+i sin /fy+i/% sin /fy+iV^} 

+«-**{~tt+*/8y}] • (20). 

To obtain the mean velocity parallel to x of a particle, we must add to (25), the 
terms previously investigated and expressed by (20). If we call the total u, we have 

„' = V™^^ . . . (27). 

At a short distance from the bottom e"^ becomes insensible, and we have simply 

^ = |^ s ^ 2 ^ e --%(l--2%) (28), 

^ = _?^^V^(--lt+f^). ...... (29). 

The steady motion expressed by (28) and (29) is of a very simple character. It 
consists of a series of vortices periodic with respect to a? in a distance |X. For a given 
x the horizontal motion is of one sign near the bottom, and of the opposite sign at a 
distance from it, the place of transition being at y==(2&)"~ 1 ==\/47r. The horizontal 
motion of the first order near the bottom being by (18) it=u Q cos hx cos nt, we see 
that it is a maximum when kx=0, ir 9 2tt, ... If we call these places loops, and the. 
places of minimum velocity nodes, (29) shows that v is negative and a maximum at 
the loops, positive and a maximum at the nodes. The fluid therefore rises from the 
bottom over the nodes and falls back again over the loops, the horizontal motion near 
the bottom being thus directed towards the nodes and from the loops. The maximum 
horizontal motion is simply f % 3 /V, and is independent of the value of v. We cannot, 
therefore, avoid considering this motion by supposing the coefficient of viscosity to be 
very small, the maintenance of the vortices becoming easier in the same proportion as 
the forces tending to produce the vortical motion diminish. 

To ascertain the character of the motion quite close to the bottom, we must include 
the terms in e~^. When y is extremely small 

u , =u sr V" 1 sin 2Tcx{— ifiy+ . . *} ....... (30), 

so that the motion is here in the opposite direction to that which prevails when e~~^ 
can be neglected. 
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A few corresponding values of fiy and of —- (sin /3y+ ih ^)<3" s5y +f are annexed, in 
order to show the distribution of velocities within the thin frictional layer. 
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It appears that (sin 2kx being positive) the velocity is negative from the plate out- 
wards until fiy somewhat exceeds \ir, after which it is positive, until reversed by the 
factor (1 — 2%). The greatest negative velocity in the layer is about \ of that which 
is found at a little distance outside the layer. 

Faraday found that fine sand, scattered over the bottom, tended to collect at the 
loops. This is in agreement with what the present calculation would lead us to 
expect, provided that we can suppose that the sand is controlled by the layer at the 
bottom whose motion is negative. The exceeding thinness of the layer, however, pre- 
sents itself as a difficulty. The subject requires further experimental investigation ; 
but in the meantime the following data may be worth notice, though in some respects, 
e.g., the shallowness of the liquid in relation to the wave-length, the circumstances 
differed materially from those assumed in the theoretical investigation. 

The liquid was water (z/=*014 O.G.S.), and the period of vibration was -f^, so that 
n=2nX 15. The thickness of the layer 



fTT* 



V2v 



'0135 centim. 



Measurements of the diameters of the particles of sand gave about '02 centim., so 
that the grains would be almost wholly immersed in the negative layer, even if isolated. 
It seems therefore that the observed motion to the loops gives rise in this case to no 
difficulty. But it is possible that the behaviour of the sand is materially influenced 
by the vertical motion of the vessel by which in these experiments the liquid vibra- 
tions are maintained.* 

§ 2. In the problem to which we now proceed the motion will be supposed to have 
its origin in the assumed motion of a flexible plate situated when in equilibrium at 
y=0. Thus for a first approximation we take i£=0, v=v sin Jcxe int , when y=0, and 
the question is to investigate the resulting motion of the fluid in contact with the 
plate. 

* See a paper " On the Crispations of Fluid resting upon a Vibrating Support," Phil. Mag., July, 1883. 
MDCCCLXXXIV. C 
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The solution to a first approximation is readily obtained. As in (10), (11), we have 

\Jt = \fi 1 +\l*2 = e^-co8Jcx(Ae-*y+Ce-'fy) (31), 



in which we may take as before 



V= A/~(l+i)=l3(l+i) (32). 



By the condition at y=0, 



A = -! c > 






so 



that 



v e int cos j cx r y I 



• ♦«#«« 



(33), 



^ 



Vr£ int COS fa f 7 , , 7 , ,, 1 , , 

= fe-y U^ y -^4 . (34). 



In passing to real quantities it will be convenient to write 



v, 



h-h 



\,=H^ 



• *> 



C V 



(35). 



Thus throwing away the imaginary parts of (33), (34), we get 



\jj= cos fee < 



— ^y-e""^cos (^+e+^7r)+e~^cos (nt+e—fiy) 



Jc 



(36), 



u=</2.pll cos Jcx< e~*y cos (^+e+|7r)— e~^cos (wt+c+iTr— y8y)| . (37), 



i?=:H sin fee 



■£v/2 6"^ cos (^+c+^7r)+ Ae""^ cos (nt + e—fiy)\ . . (38). 



From (32), (35), the approximate value of H is — v Q //5^/2 9 and that of e is -~^7r. 
More exact values will however be required later. We find 



H= 



V, 










e • r 



. (39), 



cos e: 



p-k 



=~ 1 



V{(/3-^+/3 3 } v/2 



2A 



• • 



• (40). 



The values of ^ and v above expressed give u=0 9 v=v Q sin kx cos nt, when y=Q- 
This is sufficient for a first approximation, but in proceeding further we must remember 



OF AIR OBSERVED IN KUNDT'S TUBES. 11 



that these prescribed velocities apply in strictness not to y = 0, but to y= — sin hx sin nt. 
Substituting the latter value of y in the expressions (37), and (38), we find 

w= A /2.jSH cos hx \ —hy cos (^£+ € +i^) + v^./fy cos (n£+ e+ib 7 ") f 
/3\~K . 7 . f h 



n 



sin 2hx sin ntl — ^-^-cos (^+€'+i^)+ cos (^+€+^7r) V 



y8\H . 



2n 



sin 2fee 



ft 



— -r 2 sin (e+^Tr)— sin (e+^rr) 



► + terms in 2nt. 



The first term within the bracket is of the second order in JcJ/3 relatively to the 
latter term, and may be omitted. Thus 

u= — , sm 2/cx cos €. 
2w 

The terms in 2nt we need not further examine. From (39), (40), H cos €=-— v /2/3, 
very approximately, so that we may write 

Bv 3 
w=— ^- sin 2fe . . (41). 

To the same degree of approximation, v=v sm hx cos nt, simply. 

We have next, as in the first problem, to consider the complete equation 

^ = «^.**±s (42) 

v v v> dxdt^v* dydt ' ' .' V ' 



in the right hand member of which we use the approximate values given by (36), 
(37), (38). Thus 

-—"=• — nH cos hx e~^ y sin (nt+e—fiy), 
and (42) becomes 



^ = ?Mi!^ . (43). 

It will be found presently that the term divided by h disappears from the final 
result, and thus we have to pursue the approximation further than might at first 
appear necessary. We may however neglect terms of order P//3 3 , in comparison with 

the principal term. Thus v 4 may be identified with — , and the equation becomes 

c 2 
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dty nkW sin 2kx eT™ f /2/3 



dtfyj 



4v*j3* 



h 



— 1 ) sin /3?/— cos fiy — 2/3y sin /3y + 2<?* 



-|3y 



whence 



,/,= 



«i:H 3 sin 2&c e Pv 



i. 

2k 



L f ) sin /3y+f cos $y+\$y sin /Jy+£e to 



And 



u- 



dty rikW sm2kxe. ~ Py 



(44), 



(45). 



- 2 ) ainfiy— rv cos/Sy—±Py sm(3y+±py cos py—\e-* 



M 



21c 



(46 j. 



To obtain the value of w at the surface of the plate it will be sufficient to put y=0 
in (46). Thus 



u- 



nkW sin 2fa 
4z^/3 3 



2fc 4 



(47). 



By (32), (39) 



v. 



lc 



4^ 3 /3 2 2»\ i ~,3/ 2VV """V' 



if as before we put V for &/w. Thus in (47) 



u 



=w(~k~*) ski2kx ' ' ^ 48) - 



To obtain the complete value of w at the surface of the plate, corresponding to (37), 
(46), we have to add to (48) that given in (41). The term of lowest order disappears, 
and we are left simply with 



u- 



o V 



• • • 



(49). 



In like manner we find for the complete value of v at the surface of the plate 
corresponding to (38), (45), 

llvr?k cos 27cx 



v=v () sin Jcx cos nt—~ ° n ^ Tr 

opv 



(50). 



The values of u and ^expressed in (49) and the second part of (50) must be can- 
celled by a suitable choice of the complementary function, satisfying ^%Jj=z0 9 so that 
to the second order of approximation the fluid in contact with the plate may have no 
relative motion. 

The complementary function is 



\jj=(A+By)e- u y sin 2kx, 



whence 
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v= — 2k(A+By)e- 2k v cos 2kx. 



Detennimng the constants as indicated above we get 



u=-^-(l^2ky)e-^sin2kx . . (51), 



v=—^-(ll + 6fiy)e' m ^cos2kx ...... (52), 

The velocities given by (51), (52) are the only part of the motion of the second 
order which is sensible beyond a very small distance from the vibrating plate. The 
nodes of the plate (where sand would collect) are at the points given by kx=Q, tt, 2tt . . . , 
and the loops at the points Jcx=^tt, f it . . . At the former points v is negative, and at 
the latter positive. For Jcx=^tt, u is positive, and for fe==§7r, u is negative. 



< A *~ 



V ^ 



^ A 



•I- 



-^7T 7T f 7T 

node loop node loop 

The magnitude of the vortical motion is independent of the coefficient of friction. 

The complete value of u to the second order of approximation (except the terms in 
2nt) is obtained by adding together (37), (46), and (51), and it will contain the term 
divided by h in (46), whose appearance, however, is misleading. The objectionable 
term will be got rid of, if we express the mean velocity of a particle, instead of as in 
(46), the mean velocity at a point. For this purpose we are to add to (46), (51), the 
mean value of 

gdu du 
^dm +V dy } 

as calculated from the first approximation where 

£= udt, 7j=\vdt. 

du 

As in the former problem the mean value of £— is zero. 

x dx 

du f 

Multiplying together — , and \vdt as found from (37), (38), and rejecting the terms 
in 2nt, we get with omission of £ 2 , 
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^^^^{la-ftyXsin/Sy-cos^+ie-^}. . . . (59).. 



m which we may write 

k^W k&W nkW 



n n/3 2 4v%/3 2 



Combining (53), (46), and (51), w T e get finally 



x& — — 



4^/3* 



<?' 3 sin 2&# 



j —2 sin #y— fe-* 1 +|^(l — 2ky)e" 2 ^ sin 2te 



■6"^(2 sin $/+|<T^)+f(l — 2%)e~ 2 *n (54), 



which expresses the mean particle velocity. 
When fit/ is very small, (54) gives 

A V 

from which it appears that quite close to the plate the mean velocity is in the 
opposite direction to that which is found outside the frictional layer. 

§ 3. In the third problem, relating to Kundt's tubes, the fluid must be treated as 
compressible, as the motion is supposed to be approximately in one dimension, parallel 
(say) to x. The solution to a first approximation i$ merely an adaptation to two 
dimensions of the corresponding solution for a tu^e of revolution by Kirchhoff,"* 
simplified by the neglect of the terms relating to the development and conduction 
of heat. It is probable that the solution to the second order would be practicable 
also for a tube of revolution, but for the sake of simplicity I have adhered to the 
case of two dimensions. The most important point in which the two problems are 
likely to differ can be investigated very simply, without a complete solution. 

If we suppose p=a*p, and write or for logp— logp , the fundamental equations 

are 

dor du du du , 9 , ,d (du dv\ . . 

a z ~~~=— -■- — u^r-—v ■- +^V u+v 7- - +~r (56), 

dx at dx dy dx\dx ay) v ' 

with a corresponding equation for v 9 and the equation of continuity, 

du . dv . do . do do /pH7 v 

~+T"+~r;+^~^+ / * ; ~7~ :== » . . (57). 

dx dy dt dx dy v 7 

* Pogg. Ann., t. cxxxiv,, 1868. 
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Whatever may be the actual values of u and r, we may write 

u -^±jlJ± d AJ± / 58 \ 

dx dy ay ax 

in which 

du dv „ tfw <fo , v 

'V <p="ir+"7"> V^A=i — T" (59). 

1 ax ay dy ax ' 

From (56), (57), 

2 ,d\da du 3 $w rfw ,d / <icr dcr\ //?a\ 

dtjdx dt cfo? % £&c\ <i^ <%/ * * # ^ '* 

dtj dy ^ & <% <%\ <fo? <%/ 

Again from (60), (61), 

f , d , ,d\ n d^a d f d<r . da\ , , i\ * / d<r . da s 

d / du du\ d I dv dv\ , . 



^,2? 



For the first approximation the terms of the second order in w, v, and <r are to 
be omitted. If we assume that as functions of t, all the periodic quantities are 
proportional to e in \ and write q for a?-{-inv-\-inv', (62) becomes 



Now by (57), (59), 



so 



that 



and 



%i 



gV 2 cr+?iV=0 .......... (63). 



V 2 c& = — incr =i- V 2 o\ 



T n 



n ax ay n dy ax ' 



Substituting in (60), (61), with omission of terms of the second order, we get in 
view of (63), 

(,V 3 -^=0, (vv 3 -m)j=0, 

* It is unnecessary to add a complementary function <fi\ satisfying v 2 0' = O, as the motion corre- 
sponding thereto may be regarded as covered by yjr. 
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whence 

(^V s — in)\jj—0 .......... (65). 

If we eliminate cr directly from the fundamental equations (56), we get 

, d \ , d I du , du\ d I dv , tfe\ d , , x , d , . x 

'** * t'^+j^t+jgt (66) . 



<fc c??// r dx dy 

If we now assume that as functions of x the quantities cr, \p, &c., are proportional 
to e iJcx , equations (63), (65) may be written 

/^_F 3 V=0 ? where W*=W-- . ..... (67), 

\dy J q 

^-PWo, where k'*=k*+—. ..... (68). 

• If the origin for y be in the middle between the two parallel boundaries, cr must be 
an even function of y, and xjj must be an odd function. Thus we may write 

cr=A.coahk"y.e in *cf hf 9 xjj=B sinh k'y. <? nt e itv .... (69), 



u=( — —A cosh Fy+FB sinh &'?/W e i]cx 
v=l — ~A sinh k"y—ikB sinh F^/ je^ e*** 



> 



* • » v * \ I V.' /* 



If the fixed walls are situated at y—^y^ u and v must vanish for these values of y. 
Eliminating from (70) the ratio of A to B, we get as the equation for determining k, 

W tanh kfy x = Wk" tanh k"y l . . . . . . . . (71), 

in which k\ k" are given as functions of h by (67), (68). We now introduce further 
approximations dependent upon the assumption that the direct influence of friction 
extends through a layer whose thickness is a small fraction only of y v On this suppo- 
sition Tc is large, and k" is small, so that we may put tanh k f y 1 '=- ±1, tanh Wyy=- ±&'Vi- 
Equation (71) then becomes 

K — fC rC y-i . f ........ \ / ^/j 
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or if we introduce the values of If, h" from (67), (68), 



»=[V-ty Viy /(»+? 



9 



n* %* 



Since in/v is great, P~ — =— approximately. 

q a 

Thus 



2 



and 



If we write k=k 1 -{-ik 



2> 



w 



n 



2 



Vi v ( F +7 



= 7,1 + 



a 



I 



yi 



V v j 



% 






j 



/- 



% 



a 



\f In 



"I 



*i=±r -I H — |^~j> ^ == Fr 



V 2» 






(73). 



(^), 



which agrees with the result given in § 347 (11) of my book on the Theory of Sound. 
In taking approximate forms for (70), we must distinguish which half of the 
symmetrical motion we contemplate. If we choose that for which y is negative, we 
replace cosh ley and sinh h'y by \*T Vv . For cosh k"y we may write unity, and for 
sinh k"y simply k"y. If we change the arbitrary multiplier so that the maximum 
value of u is unity, we have 

u= (— 1 + e -^y+yi)) e ^v w 

V\2fl 1 



in which, of course, u and v vanish when y-=. —y v 

If in (75) we change k into — i, and then take the mean, we obtain 



u=(—l + e-^y^) cos kx e int 

v= — - A ~+ e-^y+yJ ) sin kx e 
k'\Vi J 



> 



(76). 



Although k is not absolutely a real quantity, we may consider it to be so with 
sufficient approximation for our purpose. If we write as before 



h '= V{£)< i + i )=ft i+i )' 



MPCcclxxxiv. 
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we get from (76) in terms of real quantities 



u= cos hx[_— cos nt-\- e~^ + ^ cos {nt—^y-^y^}^ 



h 






0s/2 



sin HjQc 



Lyi 



cos (nt— ^7r)+^ {y+yi) cos {n£— ^7r— fi(y+yi)} 



> . (77). 



It will shorten the expressions with which we have to deal if we measure y from 
the wall (on the negative side) instead of as hitherto from the plane of symmetry, for 
which purpose we must write y for y-\-y v Thus 



u= cos hx{ — cos ntf+tT* cos (nt—/3y)} 






Jc sin hx 



y\-y 



Vi 



cos (nt — ^7r) — e ^ cos (nt — ■ ^77 — fiy) 



From (78) approximately 



. . (78). 



vV=i8 v /2. cos fee e"^ sin (n«— |7r— /3y) (79), 



du t dv 7.7 

— +— ==#smtecos ?i£ . 



. . . . (80), 



J^^ v J^T = Lj c j3 s i n 2fep e -^(— cos /3y+e~^)+ terms in 2n£ '. . (81), 

n • y. Cult 



dx 



dn dt) \ 

T "^~d~ r x l J= ""4^ sin 2kx e~^(sin f3y+ cos £y)+ terms in 2nt . . (82). 

As in former problems the periodic terms in 2nt will be omitted. For the non- 
periodic part of \jf of the second order, we have from (66) 

V 4 ^=~^sin2/^e^{sin/32/+3cos/32/---2e--^} .... (83). 



d* 



In this we identify v 4 with — , so that 

ay 



h sin 2hx er^J c . _ - 

^ = — 16^3 * sm ^ +3 C0S $/+i 6 ^ 



. (84), 



to which must be added a complementary function, satisfying ^xfj=0 } of the form 



sin 2to 



*fr= 16 ^ 3 ( A sinh ^^(2/i — 2/)+ B (2/i — 2/) cosh 2%!— y)} 



• • 



(85), 
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or as we may take it approximately, if y x be small compared with the wave-length X, 

+=*^{*toi-y)+w<3ri-y)*} (86). 

The value of <r to a second approximation would have to be investigated by means 
of (62). It will be composed of two parts, the first independent of t, the second a 
harmonic function of 2nt. In calculating the part of d<j>/dx independent of tfrom 

o , ' da da da 

' dt dx dy 

we shall obtain nothing from daldt. In the remaining terms on the right-hand side it 
wU. be suffice to 4% the values of u. „, „ of Z first approval. Fro m 

da du dv 

dt dx daj 

in conjunction with (80), we get 



-~ 9 sin kxsmnt, 
a 



whence 



d 2 d> hu 3 
^)i = 2^ COs3fae " %sin ^ 



It is easily seen from this that the part of u resulting from d<f>/dx is of order F/3 3 in 
comparison with the part (87) resulting from \fj, and may be omitted. 

Accordingly by (84), with introduction of the value of /3 and (in order to restore 



2 



and from (86) 



homogeneity) of u 

^= — — — {4 am fiy +2 cos fly + e-*} . ..... (87), 

2ku n $ cos 2Jcxe~P» , .«.„.«.,., , 

v=~ 8/3a {Bmpy+3coaPy+y-»} (88); 

u =- JSL ^~-'i A '+ sB '(yi--yY} (89). 

v =- 2h %T kx {^-y)+ny^-y) s } ■ ...... (90). 

When y=0, the complete values of u and v, as given by the four last equations, 
must vanish. Determining in this way the arbitrary constants A 7 and B', we get as 
the complete values at any point, 

i) 2 
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1£= 



w 3 sin 2km 
8a 



<r*(4 sin/3y+2 co$/3y+e~^)+ 



b 9.(yi-y) 9 ' 

//I 



• * i yi ij 



V: 



2fe n 9 cos 2kx ( 



j o 



(2/i—2/> 



3-1 



spa r ^ sin ^+ 3 cos £y+K*0+t£(yi-y)-"l£ 8 



l «7 Zi j< 



Outside tlie thin film of air immediately influenced by the friction we may put 
e-M=0, and then 



u _ ._ 3m q 3 sin 2fa f 1 _ c fv\ -yf 1 



16a [ 



9V 



• •*••• 



(93), 



V: 



3% 2 2& cos 2to; f (y, —yf 



16a 



* * • * • 



(94), 



From (93) we see that u changes sign as we pass from the boundary y=Q to the 
plane of symmetry y=yi, the critical value of y being y x (l — <s/i)> or *423 y x . 

The principal motion being u~-~u cos fee cos nt, the loops correspond to kx—Q, tt 9 
2tt, . . . , and the nodes correspond to \ir, f w, . . . Thus v is positive at the nodes and 
negative at the loops, vanishing of course in either case both at the wall y=Q, and at 
the plane of symmetry y=y x . 



Plane of symmetry _ 

-< — a ■ *■" - •■ •> •< — j 



Wall 



^ 



->■ 



4 ■ j "'~" ■> 



a 



-^ 



r7T 



7T 





loop node loop node 



f 71 " 



To obtain the mean velocities of the particles parallel to x y we must make an addition 
to u, as in the former problems. 

In the present case the mean value of 



du u dn 

dx dy ' 



u ( ?sin2kxe~ f * v 
ia 



e~to— cos/Jy k 



so that 



u , = _«l*^[^ {4L gin to+Se-^+i-fl 1 -*? 



8a 



Vx 



• # • (i/O )* 



hen py is small, 



u 



u * sin 2Jcx 



—-2/3?/+ • • • r (96). 



Inside the frictional layer the motion is in the same direction as just beyond ii 
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We have seen that the width of the direct current along the wall is '423 y l9 and 
that of the return current (measured up to the plane of symmetry) is '577 y l9 so that 
the direct current is distinctly narrower than the return current. This will be still 
more the case in a tube of circular section. The point under consideration depends 
only upon a complementary function analogous to (86), and is so simple that it may 
be worth while to investigate it. 

The equation for \(j is 

— — - — — 4&*f\b—0 . . . (97) 

dr 2 r dr J ™ \ '' 

but if we suppose that the radius of the tube is small in comparison with \, W may be 
omitted. The general solution is 

^={A+Br 2 +BV 3 logr+Cr 4? }sin2^ (98), 

so that 

u=\ ^={2B+B'(2logr+l) + 4Cr 2 } s i n 2 hx 9 
whence B'=0, by the condition at r=0. Again 

1 dyk 

v=z — _r. _ _ 2^{Ar~ 1 +Br+Cr 3 } cos 2kx 9 

i CMC 



whence A=0. 

We may take therefore 



u= {2B+4CV 3 } sin 2kx 
w=— • 2Jc{Br+Cr z } cos 



2hxr * * (99) * 



If v=0> when y=R, B+GR s =0, and 



w=2C(2r 3 —R 2 )sin2^ (100). 

Thus u vanishes, when 

^="72= '707 B* R-r= '293 R. 

The direct current is thus limited to an annulus of thickness *293 R, the return 
current occupying the whole interior, and having therefore a diameter of 

2X707 R=1-414R 



